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Abstract 


Recent work has shown the existence of a relativistic effect present in a single component non- 


Q. 


equilibrium fluid, corresponding to a heat flux due to an electric field LLfl. The treatment in that 
work was limited to a four-dimensional Minkowksi space-time in which the Boltzmann equation 
was treated in a special relativistic approach. The more complete framework of general relativity 
can be introduced to kinetic theory in order to describe transport processes associated to elec¬ 
tromagnetic fields. In this context the original Kaluza’s formalism is a promising approach 2-|4|. 
The present work contains a kinetic theory basis for Kaluza’s magnetohydrodynamics and gives 
a novel description for the establishment of thermodynamic forces beyond the special relativistic 
description. 
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I. INTRODUCTION 


Transport processes in fluids are describable in terms of statistical averages of micro¬ 
scopic collisional invariants which evolve according to macroscopic conservation laws. The 
mathematical form of the corresponding evolution equations is based on the knowledge of a 
distribution function and a master equation which rules its space-time behavior. Thermody¬ 
namic variables such as internal energy, heat flux and hydrodynamic velocity are identified 
with moments of the distribution and their expressions satisfy a set of partial differential 
equations that depend on the regime of validity of the distribution function. In the case of di¬ 
lute gases, the Boltzmann equation has successfully described transport processes within the 
Chapman-Enskog expansion of the distribution function in terms of powers of the Knudsen 
number Jj|. 

Nevertheless, when the gas is under the action of a field ,, the fluxes caused by the gradients 
inherent to it may be treated in the kinetic theory formalism in more than one manner. In 
the standard classical formalism, the individual particle’s acceleration is coupled to the 
fields through Newton’s law if = F e /m. In this case, the Boltzmann equation for a single 
component system reads 




( 1 ) 


Equation ([[]) is an evolution equation for the distribution function / ( x e , v e , t'j which is 
the density of gas particles in phase space. Here i = 1,2,3 and Einstein’s convention for 
equal indices has been used as in the rest of the work. The components of the external 
force are given by F e and m is the mass of the individual molecules. The term of the right 
hand side includes the changes in the number of particles in each cell of phase space due 
to collisions and its argument f f that it is a function of the product of the distribution 
functions before and after the interactions. 

If the force included in Eq. (JT|) is conservative, F e can be written as the gradient of 
a scalar field if such that v l oc iff (the comma denotes spatial derivatives). Thus the 
spatial inhomogeneities of the field constitute, in principle, thermodynamic forces capable 
of creating non-equilibrium fluxes. These thermodynamic forces are responsible of well- 
known non-equilibrium effects (Benedicks, Thomson, etc). It is interesting to notice that 
in the case of a single component system those fluxes vanish when the Chapman-Enskog 
expansion is used to solve Boltzmann’s equation in the Navier-Stokes regime. This fact will 


2 



be reviewed in section 2 of this work. 


This problem is partially overcome if the first term of Eq. (OQ) is arbitrarily suppressed 
6|, |7] or if the fluxes are established for a binary mixture in the Lorentz gas case |5|, [8]. Simple 
effects such as electric conduction and thermoelectricity require one of those approaches in 
order to describe experimental results for monocomponent charged dilute gases. 

In this work we propose the use of a different treatment of the Boltzmann equation in 
the spirit of the general theory of relativity. In this approach, individual particles “fall 
freely” in space-time so that curvature constitutes a thermodynamic force in addition to 
local variables’ gradients (VT, Vn, V«). In this type of formalism the term of Eq. 

(JIT) can be suppressed for structureless particles and the effect of “external forces” becomes 
present in the covariant derivatives that are necessary to guarantee the covariance of the 
transport system. 

Although it is well known how to describe the motion of freely falling particles in terms of 
geodesics in the case of gravitational fields = — T i ab v a v b ), the introduction of the concept 
of “free-fall” in the case of an electric held is rather subtle. One elegant way to establish the 
equation of motion of a simple particle imbedded in an electric held in terms of geodesics was 


proposed by T. Kaluza back in 1921 Q. Kaluza’s idea of a five-dimensional space-time allows 
to obtain Maxwell’s equations using Einstein’s held equations and to reproduce the Lorentz 
force exerted on a charged particle in terms of a geodesic on a curved space-time in the 
presence of source charges and currents. Decades later, the classical Kaluza framework was 
introduced in a phenomenological attempt to rewrite the equations of magnetohydrodyamics 
[3!j • The present work provides a kinetic theory support of Kaluza’s macroscopic hydrody¬ 
namics in the presence of an electric held, yielding also a novel description of old effects 
relating forces and fluxes present in the transport theory of monocomponent dilute gases. 

We have structured this paper as follows. In Sect. II we review the standard treatment of 
external forces for monocomponent dilute gases emphasizing how external forces effects are 
canceled in the hrst order in the gradients correction of the distribution function (Navier- 
Stokes regime). Also the usual argument regarding electrical conduction within kinetic 
theory in the stationary approximation is critically examined. In Sect. Ill we establish 
the system of relativistic transport equations for a gas in a hve dimensional space-time. 
In Sect. IV we establish the distribution function to hrst order in the Knudsen parameter 
following the Chapman-Enskog method and calculate of the electrical conductivity coefficient 
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in Kaluza’s formalism. Our final remarks are included in Sect. VI. 


II. THE STANDARD KINETIC THEORY APPROACH 


In this section we show how standard kinetic theory to first order in the gradients, yields 
no electrostatic contribution to the non-equilibrium distribution function for the single com¬ 
ponent charged gas. The non-relativistic Boltzmann equation for such a system in the 
presence of an electrostatic potential </> is given by 


df , _ as , df /-/(»> 

pr + v • H-V0 • — —- 

ot or m ov r 


( 2 ) 


where a relaxation time approximation has been introduced on the right hand side. Here 
f 1 ' 0 ' 1 is the local equilibrium solution, v the molecular velocity, r a relaxation time and m 
and e the mass and electric charge of the molecules respectively. For the first order in the 
gradients deviation from equilibrium / (1 ^ = / — /(°) one obtains 


f (D (df {0) , - <9/ (0) e dfW\ 

; = - T (pr +,, '^F + m v ' # "^p)- 


( 3 ) 


where /d) is given by Maxwell-Boltzmann’s distribution function. The electric field contri¬ 
bution to /d) arises from two terms, namely the last term on the right hand side of Eq. ([2]) 
and the one that appears in the time derivative since 

df d) 9/d) Qn 9/d) 9T 9/d) du 

dt dn dt dT dt du dt 

According to the Chapman-Enskog method, ^ can be expressed in terms of spatial gradients 
and external forces in the Euler regime. Calculating explicitly both contributions one obtains 

9 /(°) 


-V0 

m 


a, =-iV l0,v ^’ 


and 




(5) 


( 6 ) 


where k = v — u is the usual peculiar velocity, u being the hydrodynamic (bulk) velocity 
of the fluid. From Eqs. (J5J) and (J6]) it is clear that this standard procedure yields no 
electrostatic contribution to fd\ the reason being that the drag due to the electric force on 
the individual particles and on the fluid as a whole cancel each other out. This does not 
occur in the mixture case. 
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Alternatives to obtain a dissipative electrical current for the simple charged fluid within 
kinetic theory have been put forward, for example in Refs, jcl] and [ 7 ]. Both authors consider 
a steady state case by neglecting the time derivative term in Eq. ([2]). However it is important 
to notice that once the Chapman-Enskog hypothesis is introduced, this term is of the same 
order in the Knudsen parameter as the rest of the terms in Eq. ([2]). Indeed, Eqs. (H)) to (121) 
show that this term is of first order in the gradients as the other two terms on the right hand 
side of Eq. (J2]) are. In the following sections of this work, a general relativistic formalism is 
developed in order to establish a constitutive equation for the electrical current using kinetic 
theory to first order in the gradients without introducing approximations in the Boltzmann 
equations foreign to the Chapman-Enskog scheme. 


III. TRANSPORT EQUATIONS IN KALUZA’S FRAMEWORK 

Kaluza’s original work establishes the possibility of including electromagnetic effects 
in the geometry of space-time by showing that a charged particle moving along a geodesic 
in the proposed five-dimensional metric follows the same trajectory as the one dictated by 
the Lorentz force. Such a condition is satisfied, to first order in the components of the 
electromagnetic four potential, if the metric is given by 

ds 2 = rj^dx^'dx 1 ' + C.A^dx^dx 5 + (dx 5 ^j , (7) 


where /i and v run from 1 to 4 and ris the cartesian Minkowski metric tensor. The 
additional elements of the metric tensor are given by a normalized electromagnetic potential 
where 7U = A, 6/c , with A and 6 being the usual vector and scalar potential 

Q. 


and £ = 


= ,/ 167i f e Q 


In this scenario, the non-vanishing Christoffel symbols are 


p/' _ pM _ ip/t 

v 5 2 " ’ 


■e 


( 8 ) 


and 

1= -fi?, (9) 

where F^ w = A /i; „ — A v . jfl = A^ u — A viAt is Faraday’s field tensor. 

The equation of motion for a particle in the space-time given by Eq. © is determined 


by the geodesic equation 


dv^ 

dr 




aX 


v a v x 


0, 


( 10 ) 
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where v 11 = and r is the particle’s proper time. Comparison of Eq. (jUJ) with the equation 


dr 


of motion for a charged particle in the presence of an electromagnetic held F^ v 

e 


dr 


m 


-v a F» a , 


( 11 ) 


gives rise to the relation proposed by Kaluza: v 5 = Also, the so-called cylindrical 
condition that all derivatives with respect to x 5 vanish is imposed in the original formalism. 
This condition leads to the concept of periodic, compact extra dimensions [9]. 


For the kinetic theory purposes, the charged fluid’s state variables are given as 

W = ~ f i)fd*v, 

n . 


( 12 ) 


where 0 is taken as nv M for the particle four-flux {n being the number density) and as mv^v 1 ' 
for the energy-momentum tensor, for details on the nature of these quantities the reader can 


refer to Ref. 


lOj. The average is performed over the three dimensional velocity space, the 


corresponding invariant volume element is given, as in the special relativistic case, by 

j* 5 d W 4 ,3 

dv = 7 = l{w) d w . 


(13) 


where w l is the vector three-velocity and 7 ^) = (l — ^ the usual Lorentz factor. A 

simple derivation of this volume element can be found in Appendix B of Ref. Jll|. 

Balance equations for state variables can be obtained in this framework by multiplying 
the relativistic Boltzmann’s equation 

df 


= j un 


(14) 


by a function of the molecular four-velocity -0 (u M ), and integrating in velocity space. Here 
a comma denotes a partial derivative and the dot corresponds to the total derivative 

/ dA v x 


yR = = v» 

Jr* 


+ r 


(15) 


ydx» ' , 

where a semicolon is used to indicate a covariant derivative. The distribution function / cor¬ 
responds to the particle density in three dimensional phase space such that / vd 3 xd 3 v 
yields the number of particles in a volume d 3 xd 3 v. The collisional kernel on the right hand 
side of Eq. (1T41) accounts for particles entering and leaving a cell in phase space due to 
molecular interactions. For the second term in the left hand side of Eq. (1T41) one notices 
that 


yfj. — y a v ^ a — y a 


dv M 
dx a 


+ r »v x 


(16) 
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or 


*'■=— + r>\ (17) 

which vanishes since particles move following the geodesic equation given by Eq. (0|). This 
leads to the force-free Boltzmann equation 


^ = J <T), 


(18) 


which is consistent with the present approach where the effects of fields must be included in 
the geometry of space-time. 

Multiplying Eq. (fT8l) by a collisional invariant and integrating in velocity space leads to 


fif(v a )v fJ 'd*v 


= 0, 




where we have used the relation 


(19) 


K >= KV>(l>“)];„/ + /,«• 


( 20 ) 


In Eq. (1201) the first term on the right hand side vanishes since if only depends explicitly on 
the molecular velocity. This procedure leads to the usual balance equations for the conserved 
fluxes: 

JVJ = o, T“" = 0, (21) 

with AT and T afl being the particle four-flux and energy-momentum tensor respectively, 
given by 

AT = J fv fl d*v, (22) 

and 

= m J fv»v a d*v. (23) 

In particular, in a local equilibrium situation, the conserved fluxes are 


AT = nW, 


(24) 


and 


V 1U = pWU v +pg^. 


(25) 


where U u is the hydrodynamic four velocity and p = (ne + p) /c 2 with £ and p being the 
internal energy and hydrostatic pressure respectively [l(3|. In such regime Eqs. (12TT) evaluated 
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for / = /(°) correspond to Euler’s equations for a charged single component fluid in the 
presence of an electrostatic field. This set can be written as 


U a n a = -nU £, 

(26) 


(27) 

/ WU U \ 

(28) 

m v K = - (jr + — j p,„, 


where we have introduced the specific heat at constant number density defined as C n = 
• In Eq. (l28p . which corresponds to the momentum balance equation for /i = 1,2,3, 
the electric force acting on the fluid emerges as a consequence of curvature through the 
Christoffel symbols included on the left hand side where 

f)i 

U v Uf v = U v — + Tf x U IJ U x . (29) 

Equations (I26H28P are thus the set of transport equations for the charged single fluid in 
Kaluza’s framework in the local equilibrium, Euler regime. The momentum balance given 
by Eq. (1281) will be used in the next section in order to calculate the contribution of the 
electrostatic field to the distribution function and thus to the dissipative effects. 

IV. THE DISTRIBUTION FUNCTION IN MARLE’S APPROXIMATION 

In this section, the Chapman-Enskog method of solution for the Boltzmann equation is 
used in order to obtain the basic structure of the out of equilibrium distribution function to 
first order in the gradients, including a term dependent on the electrostatic field. In order 
to accomplish this task, we start by stating the Chapman-Enskog hypothesis, namely 

/ — / (0) + e / (1) + • • •, (30) 

where e is Knudsen’s parameter which measures the degree in which the distribution function 
deviates from equilibrium accounting for the relative scale in which macroscopic gradients 
are significant to the microscopic length. Thus, the first order in e approximation in Eq. (1301) 
corresponds to the first order in the gradients scenario, namely the Navier-Stokes regime. 





Here the equilibrium distribution function /(°) is given by the relativistic Maxwellian (or 
Juttner function) namely, 


n 


-e zc 


(31) 


Attc 3 zK 2 ( 7 ) 

where f3 = 1,2, 3, 4. In Eq. (|3Th K 2 is the modified Bessel function of the second kind and 
z = ^7 is the relativistic parameter that compares the thermal energy with the rest energy 
of a single particle of mass m. 

Moreover, in order to establish the structure of in terms of the gradients, a relaxation 
time approximation is appropriate since the flux-force relations in the case of conservative 
forces is retained in such model. The complete Kernel would only add precision to the 
transport coefficients which is beyond the task at hand. Thus, for the sake of clarity and 
simplicity, we consider Marie’s approximation and write 

/ - / (0) 


Jiff ) = 


r 


(32) 


where r is the relaxation parameter which scales as the microscopic characteristic time 
Introducing Eq. (l24ll in Boltzmann’s equation with the approximation given by Eq. 
leads to 


;f m = 


-TV r 


0f(°) df w 
J n u + 4—T u + ^—U c , 


(33) 


dn ' dT ' du a 
where the functional hypothesis, stating that the distribution function depends on space and 
time only through the state variables, has been also used. The derivatives of the distribution 
function can be readily obtained and are given by 

0) f(0) 


dn 


n 


(34) 


Qf(0) ~ y(0) (^ + U P Vf} _ g 


dT 


zc“ 


and 


df^ v, 


u a f ( 0 ) 
9 J ) 


which, when introduced in Eq. 

e/W = -r^/ (0) 


dU a zc 2 
lead to 


1 _ , (, V-Uv S(l 
n n ’“ + 1 “ 17 7 


(35) 


(36) 


\t u + 

zc A 




(37) 
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where Q (±) = K 3 (A) /K 2 (i). 

Dissipative fluxes arise from the out of equilibrium part of the distribution function fY) 
and are coupled linearly, in the Navier-Stokes regime, with the thermodynamic forces of 
the same tensorial rank. These forces are given by the gradients of the state variables. In 
particular, the Chapman-Enskog procedure demands, for existence of the solution, imposing 
the lower order balance equations as a way of writing the state variables’ time derivatives, 
contained in the four-gradients in Eq. (DTI) , in terms of the corresponding spatial gradients. 
Thus, in the Navier-Stokes regime, the Euler local equilibrium equations Eqs. (12011281) . need 
to be introduced in Eq. (J37J). In order to perform such a task it is convenient to separate 
time and space derivatives in each term. Also, since e/^ is an invariant, it can be calculated 
in the fluid’s comoving frame. The corresponding equation, which is derived in detail in the 
Appendix, reads 


e/ 


(i) 


-T7 (fc)/ (0) 






h 

T 


+ 


7 (fc)fco^ 
zc 2 



(38) 


where £ = 1,2, 3. Here k £ is the chaotic or peculiar velocity which corresponds to the velocity 
of the molecules measured in the comoving frame. 

The last term in Eq. ()38l) includes the contribution of the electromagnetic potential to 
the distribution function, as emphasized in the previous section. Focusing only on that term 
one obtains (£ = 1, 2, 3) 

07/o f)lJ a 

k a k»U« = Kk l -^j + k a — + r^W'Wka. (39) 

Notice that this expression contains two contributions. The last term on the right hand side 
of Eq. (1391) is due to the effect of curvature on the motion of the molecules. Meanwhile, the 
second term is also affected by the electromagnetic field since 

k a ^ ~ (g av + U ~^) kaP,u ~ T a uX U x U^k a , 

where Eq. (1281) has been introduced. Notice that this effect corresponds to the electromag¬ 
netic field acting, through the curvature of space-time, on the motion of the fluid as a bulk. 
These two terms, both proportional to the components of the Faraday tensor, are present in 
ef^ and are thus capable of giving rise to dissipation. In the case of a purely electrostatic 
field and since all terms are calculated in the comoving frame one obtains 


Y a vX U x U v k a = 2 


qc 


r Uke, 


(40) 
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and 


T° x U x k u k a = 0. (41) 

Introducing these results in Eq. (1381) and defining as ef^ the part of the distribution 
function depending on the electric field one obtains 


which, using Eq. (J8|) yields 


e/ " = 1“ 7 (U ((,,I 1A. 


e/t 0 = T -^l { °hlk) k lV- 


(42) 


(43) 


Equation (1431) shows that in Kaluza’s framework the non-equilibrium part of the distribution 
function features a contribution arising from the electric field which leads to dissipative fluxes 
in the single component fluid. As mentioned in Sect. I, the establishment of an electrical 
conductivity linked to the electrostatic potential through Ohm’s laws is of particular interest. 

The general expression for the electrical current in relativistic kinetic theory is given by 
the average m 


jP = f (r/(k)k p ) ef^d'K. 
The integral can be readily evaluated using 


(44) 


d* K = 47rc 3 y / 7 ( 2 fc) - 1 dj, 


and 


■„2 _ „2 | ^(fe) 


k l = c" 


2 -1 


T?*) 


(45) 


(46) 


Using the previous equations one obtains the following expression for the electrical current 


J‘ = r^g C- 

m \z 


The heat flux in Eq. (1471) can be written in the low z limit as 


(47) 


r ~ t 


nq 


m 


1 + y+ 


(48) 


which is consistent with the non-relativistic results obtained in Ref. 6, 7]. 
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V. FINAL REMARKS 


The concept of external forces is foreign to the spirit of the theory of general relativity, 


external forces are replaced by space-time curvature. Following this idea, it is natural to ex¬ 


plore the possibility of introducing space-time curvature into kinetic theory suppressing the 
external force term acting on single particles in Boltzmann’s equation. While the introduc¬ 
tion of curvature in the gravitational held case is rather straightforward, the electromagnetic 
counterpart has been an open question since the introduction of general relativity. One ap¬ 


proach that has proved useful in describing the dynamics of a charged particle in space time 


in terms of geodesics is the Kaluza’s original formalism first introduced in 1921. 

In this work we show how the old problem regarding how to obtain Ohm’s law for a single 


component charged dilute gas from kinetic theory consistently with the Chapman-Enskog 


expansion is solved if Kaluza’s description is applied to Boltzmann’s equation. Moreover, the 
establishment of vector fluxes such as heat and the study of thermoelectric effects through 
these techniques is feasible and will be addressed in the near future. 

In a recent publication two of us calculated the counterpart of Eq. (1T21) in ecial relativity 



introducing the gradient of the electrostatic potential as an external force |l|. In that work, 


the deviation from equilibrium due to this force vanishes in the non-relativistic limit. In 
that case, the electrostatic contribution to lowest order in 2 cancels due to the same reasons 
as the ones exhibited in Sec. II, that is the equivalency of the force acting on the fluid 
as a whole and on the individual molecules. On the other hand, as here shown, Ohm’s 
law can be obtained for a single species charged fluid identifying the electrostatic force as 
curvature. In this case the molecules move on geodesics and thus the acceleration term in 
Boltzmann’s equation is not present, leaving the geometry to act on the distribution function 
solely through the momentum balance equation for the fluid. 

As a final comment we wish to point out that the fact that a problem related with the 
electrical conductivity for a single component non-relativistic gas has been solved here using 
a relativistic formalism is somehow natural since all electromagnetic phenomena possess a 
relativistic character. Other relativistic and magnetic effects will also be accounted for in 
future work. 
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Appendix 


In this appendix, the steps leading to Eq. (1381) are shown to some detail. After Chapman- 
Enskog’s hypothesis has been introduced in Boltzmann’s equation together with the deriva¬ 
tives of Juttner’s distribution function, one obtains Eq. (1H71) in which time and space deriva¬ 
tives need to be separated. Chapman-Enskog’s method requires the introduction of Euler’s 
equations in order to write time derivatives in terms of spatial gradients in the Navier-Stokes 
regime. 

For the sake of clarity, such separation is performed only for the density and temperature 
derivatives as a first step in the main paper. In this way one can isolate the term correspond¬ 
ing to the derivative of the hydrodynamic velocity from which the electric contributions will 
arise. Also, since the calculation is carried out in the comoving frame one has U a = 0, c 


and v a = 


7(fc) 


peculiar velocity 


k, c _, where k is the velocity measured in such frame, that is the chaotic or 


lOj. Proceeding in such a way one can rewrite Eq. (1371) as follows 

v v U„ 



( 1 ) 

\ 1 dT 

2) 

J T dt 


7 lk)^k c 


zc A 


K 


(49) 


where now i = 1, 2, 3 and the greek indices /i and v still run up to 4. 

As mentioned above and in the main text, time derivatives need to be written in terms 
of thermodynamic forces by means of Euler’s equations which can be readily obtained from 
Boltzmann’s equation. Indeed, for the density one uses the continuity equation, Eq. 
which in the comoving frame yields 


dn 

~dt 


= 0 


For the temperature one uses the internal energy balance together with e >v = C n T u which 
yields Eq. (1271) . Thus in the comoving frame one also has that 

d T 


dt 


= 0 


Because of this, the first square bracket in Eq. (FUJI) vanishes so that 


£ / (1) = r / (0) 7(fc) { k 


n 


n.£ + I 1 


v v U v 


ZC“ 


5 ( 1 )' 


T 


7 {k)k^k c 


U. 


zc* 




(50) 
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which is precisely Eq. (1381) in the main text introducing v v U. v = — 7 (fc)C 2 . The subsequent 
separation and careful treatment of the last term, which involves the covariant derivative, 
is carried out to detail in Eqs. 039H43p . 
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